Clebsch-Gordan Construction of Lattice Interpolating Fields for Excited Baryons 
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Large sets of baryon interpolating field operators are developed for use in lattice QCD studies 
of baryons with zero momentum. Operators are classified according to the double-valued irre- 
ducible representations of the octahedral group. At first, three-quark smeared, local operators are 
constructed for each isospin and strangeness and they are classified according to their symmetry 
with respect to exchange of Dirac indices. Nonlocal baryon operators are formulated in a second 
step as direct products of the spinor structures of smeared, local operators together with gauge- 
covariant lattice displacements of one or more of the smeared quark fields. Linear combinations of 
direct products of spinorial and spatial irreducible representations are then formed with appropriate 
Clebsch-Gordan coefficients of the octahedral group. The construction attempts to maintain max- 
imal overlap with the continuum SU(2) group in order to provide a physically interpretable basis. 
Nonlocal operators provide direct couplings to states that have nonzero orbital angular momentum. 
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I. INTRODUCTION 

The theoretical determination of the spectrum of 
baryon resonances from QCD continues to be an im- 
portant goal. Lattice QCD calculations have succeeded 
in part to meet this goal by providing results for the 
lowest-mass baryon of each isospin in the quenched ap- 
proximation using overly large masses for the u and d 
quarks. [1, 2] 

Most lattice simulations to date have used restricted 
sets of operators appropriate for J p = l/2 ± , 3/2 ± states. 
Masses of low-lying, positive-parity baryons are repro- 
duced with approximately 10% deviations from exper- 
imental values using the quenched approximation [1]. 
Much less is known about higher-spin excited states. The 
first preliminary lattice calculation of 5/2 ± N* masses 
was reported by the Lattice Hadron Physics Collabora- 
tion [3] using one of the operators that we develop in 
this paper [4]. Results for excited baryons were also re- 
ported based on the use of different radial smearings of 
the quarks in Refs. [5, 6]. Recently, studies of negative- 
parity baryons have been reported by several groups [7- 
13]. Nemoto et al. and Melnitchouk et at consid- 
ered the A(1405) baryon, which is the lightest negative- 
parity baryon despite its nonzero strangeness. They used 
a three-quark interpolating field operator motivated by 
the spin-flavor SU (6) quark model and concluded that 
A(1405) was not evident in their lattice calculations. 

To improve upon our understanding of the resonance 
spectrum, correlation matrices will be needed, necessitat- 
ing the construction of sufficiently large sets of baryon 
and multi-hadron operators. The correlation-matrix 
method [14, 15] has been used to determine the spectrum 



of glueball masses by Morningstar and Peardon [16]. A 
large number of interpolating field operators was used 
to form matrices of lattice correlation functions. The 
spectrum of effective masses was obtained by diagonaliz- 
ing the matrices of correlation functions to isolate mass 
eigenstates for each symmetry channel. In effect, one al- 
lows the dynamics to determine the optimal linear com- 
bination of operators that couple to each mass state. A 
similar program for baryons is being undertaken by the 
Lattice Hadron Physics Collaboration. The first step is 
to determine a large number of suitable baryon inter- 
polating field operators that correspond to states of zero 
momentum, definite parity and values of angular momen- 
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Due to the complexity of the operator construction 
and the importance of providing checks on the final re- 
sults, we have been pursuing two different approaches. 
An alternative method, based on a computational imple- 
mentation of the group projection operation, is presented 
elsewhere [17, 18]. 

On a cubic lattice, the continuum rotational symme- 
try is broken to the finite octahedral group, O. States 
of definite angular momentum correspond to states that 
occur in certain patterns distributed over the irreducible 
representations (IRs) of O. Although mass calculations 
are insensitive to the spin projection J z , other applica- 
tions do require baryons with a definite spin projection. 
We develop operators that are IRs of O using a basis 
that corresponds as closely as possible to the continuum 
J, J z states in order to have operators for applications 
that require spin projection. 

It is important to use smearing of the quark fields and 
to have nonlocal baryon operators as well as the usual 
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local operators. Smeared and nonlocal operators provide 
a variety of radial and angular distributions of quarks 
within a baryon so as to couple efficiently to excited 
states. Nonlocal operators are needed in order to realize 
spins J > | and simply to enlarge the sets of operators. 

In this paper, we first review some basic facts about 
the octahedral group for integer and half-integer spins in 
Section II [19-24] and review a useful notation for Dirac 
indices based on p-spin. Two basic types of three-quark 
operators are considered: quasi-local and nonlocal. Each 
quark field in a quasi-local baryon operator is smeared 
about a common point x using the same cubically sym- 
metric form of smearing. Quasi-local operators include 
local operators as a special case, i.e., when the smear- 
ing is omitted. We develop IRs for quasi-local operators 
in Section III for each baryon: N, A, A, E, S and O. 
This amounts to determining all allowed combinations 
of Dirac indices for each flavor symmetry and classifying 
them into IRs of the octahedral group. 

The quasi-local operators provide templates that are 
used for the construction of nonlocal operators in Sec- 
tion IV. Nonlocal operators are formed by applying ex- 
tra lattice displacements to one (or more) of the smeared 
quark fields, thus providing a smearing distribution that 
differs from that used for the other quarks. In the sim- 
plest case, the combination of extra displacements used 
is cubically symmetric and only changes the radial distri- 
bution of the smearing. In other cases the combinations 
of the extra displacements transform as IRs of the octa- 
hedral group, and are chosen to correspond as closely as 
possible to spherical harmonics. The IRs of lattice dis- 
placements and IRs of Dirac indices are combined to form 
overall IRs for the nonlocal operators using an appropri- 
ate set of Clebsch-Gordan coefficients of the octahedral 
group. Some concluding remarks are presented in Sec- 
tion V. 



II. OCTAHEDRAL GROUP AND LATTICE 
OPERATORS 

In lattice QCD, hadron field operators are composed 
of quark and gluon fields on a spatially-isotropic cubic 
lattice. The lattice is symmetric with respect to a re- 
stricted set of rotations about spatial axes that form the 
octahedral group, O, which is a subgroup of the con- 
tinuum rotational group SU(2). The octahedral group 
consists of 24 group elements, each corresponding to a 
discrete rotation that leaves invariant a cube, or an oc- 
tahedron embedded within the cube. When the objects 
that are rotated involve half-integer values of the angu- 
lar momentum, the number of group elements doubles to 
extend the range of rotational angles from 27r to 47r, form- 
ing the double-valued representations of the octahedral 
group, referred to as O d . 

Spatial inversion commutes with all rotations and to- 
gether with the identity forms a two-element point group. 
Taking inversion together with the finite rotational group 



simply doubles the number of group elements, giving the 
group for half-integer spins. 

Given a lattice interpolating operator for a baryon, one 
may generate other operators by applying the elements 
of 0® to the given operator. This produces a set of op- 
erators that transform amongst themselves with respect 
to , and thus these operators Oi form the basis of a 
representation of the group. When a group element G a 
is applied to operator Oi in the set, the result is a linear 
combination of other operators in the set, YV OjTji(G a ), 
where Tji(G a ) is a matrix representation of the octahe- 
dral group. Such matrix representations are in general 
reducible. In order to identify operators that correspond 
to baryons with specific lattice symmetries, it is neces- 
sary to block-diagonalize Tji, each block corresponding 
to an irreducible representation of the octahedral group. 
This task is facilitated by a judicious choice of IR ba- 
sis vectors for the octahedral group, such as the "cubic 
harmonics" or "lattice harmonics" of Refs. [25, 26]. 



A. Integer angular momentum : O 

The octahedral group O has five IRs, namely 
Ai, A%, E,T\ and with dimension 1, 1, 2, 3 and 3, re- 
spectively, where we follow the conventions of Ref. [19]. 
The patterns of IRs of O that correspond to IRs of the 
continuum rotational group SU(2) with spin J are shown 
in Table I. A J =0 state must show up in the A\ IR, 

TABLE I: The subduction of SU(2) to IR A of O for integer 
J. 
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but in no other IR of O. A J = 1 state must show up in 
the T\ IR but in no other IR. A J = 2 state must show 
up in the E and T 2 IRs. 

Lattice displacements form representations of O cor- 
responding to integer angular momenta. We choose the 
standard "lattice harmonics" that are shown in Table II 
as the appropriate basis vectors because they have a 
straightforward connection to IRs of the rotation group 
SU(2) in the continuum limit. For example, the Y\. m 
spherical harmonics for m = 1,0,-1 provide a basis for 
the three-dimensional T\ IR. The same basis convention 
for T\ appears in Ref. [27], and the same basis convention 
for E appears in Ref. [28] . 

Any quantities that transform in the same fashion as 
the basis vectors provide a suitable IR for the octahedral 
group. We will show in Section IV how to use combina- 
tions of lattice displacements of quark fields in order to 
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TABLE II: Basis of irreducible representations of in terms 
of spherical harmonics, Y l>m for the lowest values of I. is 
the dimension of the IR. The lattice harmonics are understood 
to be evaluated on a cubic lattice. 

A d,A row 1 row 2 row 3 

Ai 1 Yo,o 

A 2 1 ^5(^3,2-^3,-2) 

E 2 y 2 , ^(Y2,2+Y2,-2) 

Ti 3 Yi,i Fi.o Fi,-i 

t 2 3 y 2 ,i ^5(^2,2-^2,-2) y 2 ,-i 



realize the same transformations as the "lattice harmon- 
ics". 



B. Half-integer angular momenta: O 

The eight IRs of the double-valued representations of 
the octahedral group, O d , include A 1 ,A 2 ,E,T l ,T 2 for 
integer spins and Gx,G 2 , and H for half- integer spins. 
The additional IRs Gi,G2 and H have dimensions 2, 2, 
and 4, respectively; these are the appropriate IRs for 
baryon operators on a cubic lattice. Table III shows the 
patterns within O d that correspond to some half-integer 
values of J. For example, a J =1/2 baryon state should 

TABLE III: The subduction of SU(2) to IR A of O d for half- 
integer J. 
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A 


1/2 


Gi 


3/2 
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5/2 


H®G 2 


7/2 


Gi®H®G 2 



show up in IR Gi. A spin 3/2 baryon should show up in 
IR H. A spin 5/2 state should show up in IRs H and G 2 
but not in G\. A J =7/2 state should show up once in 
IRs d, H and G 2 . 

A suitable set of IR basis vectors for half-integer an- 
gular momenta is given by the eigenstates | J, m) of J 2 
and J z that are listed in Table IV. Explicit forms of the 
Gi and H basis states for products of three Dirac spinors 
are given in Appendix B. Note that the G 2 basis cannot 
be built using three Dirac spinors unless orbital angular 
momentum is added. 



C. Smearing and smearing parity 

The first step in the construction of field operators suit- 
able for baryons is to specify primitive three-quark op- 



TABLE IV: Correspondence of our choice of rows in the Gi , 
G2 and H IRs to the eigenstates |J, m) of J 2 and J z . 
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erators. Consider a generic operator formed from three- 
quark fields as follows, 

e Qbc ^(x,^{ 2 (x,^{ 3 (x,t), (1) 

where a, b and c are color indices, fx, f 2 and fs are flavor 
indices and /U2 and ^3 are Dirac indices with values 
1 to 4. The operator is antisymmetrized in color by the 
e a b c factor when the (implicit) sums over a, b and c are 
performed. 

The use of gauge-covariant quark-field smearing, such 
as Gaussian smearing [29], Jacobi smearing [30] or so- 
called Wuppcrtal smearing [31], is important for en- 
hancing the coupling to the low- lying states. Gauge- 
link smearing [32-34] further reduces the coupling to the 
short-wavelength modes of the theory. Schematically, the 
smearing replaces each unsmeared field by a sum of fields 
with a distribution function as follows, 

g M (x,t) = 5^G(x,x + y)«„(x + y,t), (2) 

y 

where (x, t) denotes an unsmeared field at point x and 
the smearing distribution function G is gauge covariant. 
When the smearing distribution is cubically symmetric 
about point x and is the same for each quark field, the 
baryon operator of Eq. (1) is referred to as quasi- local. 
Quasi-local operators have the same transformations un- 
der the octahedral group as unsmeared operators. 

Nonlocal operators differ because the smearing distri- 
bution of one or more quark fields is altered by extra lat- 
tice displacements. An example is the covariant deriva- 
tive formed by a linear combination of two displacements 
of a smeared quark field, 

A^(x,i) = f/i(x)^(x + M)-C//(x-i)^(x-i,t), (3) 

where the color indices are suppressed. Equation (3) de- 
fines a new smearing distribution that is odd with respect 
to an inversion about point x. Thus, smearing can con- 
tribute in a nontrivial way to the behavior of the field 
with respect to inversion. This we call smearing parity. 

D. Inversion, Parity and p-parity 

The improper point groups Oh and consist of ro- 
tations that leave the cube invariant together with the 
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spatial inversion. The parity transformation of a Dirac 
field involves multiplication by the 74 Dirac matrix in 
addition to spatial inversion as follows, 

7>«(x,t)7>- 1 =74«(-x,t), (4) 

where V is the parity operator. Throughout this work we 
employ the Dirac-Pauli representation for Dirac 7 matri- 
ces for which 74 = diag[l, 1, — 1, — 1]. However, our re- 
sults may be used with any representation of the Dirac 
7 matrices by applying the appropriate unitary transfor- 
mation as discussed in Appendix A. 

It is convenient to express the Dirac matrices as a di- 
rect product of the form SU(2) p <g> SU(2) S where the 
SU(2) components are generated by the 2x2 Pauli ma- 
trices for spin s and p-spin p [35, 36]. See Appendix A 
for details of the construction. 

Expressed in terms of the SU(2) P (g> SU(2) S matrices, 
74 = P3 <S> 04 where p 3 = diag[l, —1] and 174 = diag[l, 1]. 
Similarly, the Dirac matrix er 2 i = §[72,71] = P4 ® 03 
where pi = diag[l, 1] and er 3 = diag[l,— 1]. With these 
conventions, a fermion field satisfies 

74<? At (x,i) = (p 3 ® (T4)? M (x,t) = p%(*,t), (5) 

and 

CT 2 iq M (x,t) = (p 4 <g> CT 3 )^(x,t) = s<^(x,i), (6) 

where Table V provides the p and s values. Thus, 
the Dirac index p = 1,2,3,4 is equivalent to a two- 
dimensional superscript corresponding to p-spin (p = 
+1, -1) and a two-dimensional subscript corresponding 
to spin (s = +1,-1), and the field may be written as 
q§(x.,i). We refer to the p value as p-parity because of 
its role in the parity transformation of Eqs. (4) and (5). 

TABLE V: Translation of the Dirac index p to p- and s-spin 
indices. Index p is expressed in the Dirac-Pauli representa- 
tion. 



Dirac index p 


P 


s 


1 


+ 


+ 


2 


+ 




3 




+ 


4 







The parity transformation of a smeared quark field can 
differ from that of an unsmeared field because the smear- 
ing parity enters. This is most easily seen by using free 
fields for which the gauge link variables are unity. Then 
the smearing distribution does not depend on the point 
x and reduces to the set of coefficients c(y) that weight 
the fields at points y away from the central point x, i.e., 
the smeared quark field is 

g M (x,t) = 5^c(y)g M (x + y,t) (7) 
y 



and the smearing parity is defined by 

c(-y)=jw(y), (8) 

where p = +1 or — 1 for even or odd smearing parity, 
respectively. The parity transformation of a smeared field 
is 

Vq^t)V~ x - ]T C (y)74^(-x-y,t) 

v 

= J4j2c(-y)q p {-x + y,t), 
y 

= 74P^ c ( y )'^(~ x + y ' i )' 

y 

= 74P«a»(- x .*)i ( 9 ) 

where the second line involves the relabeling y — ► — y 
and the third line uses the symmetry of the smearing 
distribution under inversion of y. When the gauge links 
are included so as to obtain a gauge covariant smearing a 
similar result is obtained, which holds as an average over 
gauge configurations. 

The parity transformation of a product of three 
smeared quark fields is 

P1P2P3 P1P2P3 < (-x, t)qH (-X, t)qP/ 3 (-X, t), (10) 

where we have used the notation q§ (x, t) in place of 
qfj,(x,t) and evaluated the 74 matrices using Eq. (5) to 
obtain the product of the three p-parities. The product 
P1P2P3 is referred to simply as the p-parity of the opera- 
tor and the product P1P2P3 is referred to as the smearing 
parity of the operator. 

The field operator at an arbitrary point x does not 
have a definite parity. However, in correlation functions 
projected to zero total momentum, the x dependence is 
removed by a translation following insertion of a complete 
set of intermediate states, e.g., 

C(t) = < 0|B(x,t)B(0,0)|0 > 

X 

= < 0|B(0,0)|n > e- M "* < n|B(0,0)|0 > . 

n 

(11) 

Thus the zero-momentum correlation function has 
baryon operators only at point x = where the operator 
has parity given by the product of p-parity and smearing 
parity, i.e., P1P2P3 P1P2P3 in Eq. (10). The parity of in- 
termediate state n must be the same in order to have a 
nonvanishing coupling. 

Rotations of a quark field are generated by the Dirac 
matrices ct^ = —p^®ak where indices i, j and k are cyclic 
and take the values 1, 2 and 3. Rotations are diagonal in 
p-spin and thus give a linear combination of fields with 
different s labels but the same p-parity, 

U{R)qP(^t)uHR) = J2t'(R^t)Ts's(R), (12) 

a' 



■5 



where T s i s (R) is a representation matrix of rotation R. 
This insight into the transformations of Dirac indices 
with respect to rotations is the first reason that we find 
the p, s labels useful. 

Note that a "barred" field transforms in the same way 
as a quantum "ket" when the unitary quantum operator 
U (R) is applied, i.e., 

U{R) \a) = I s ') ( s '\ U ( R ) I s ) = E I s ') T '''(R)- (13) 

s' s' 

However, "unbarred" fields also are required. Although 
they are independent fields in the Euclidean theory, their 
transformations are similar to those of quantum "bra" 
states, 

U(R)q^t)uHR) = J2q' s ,(R- 1 x,t)Tl(R). (14) 

s' 

In this paper we state results generally in terms of 
"barred" fields in order to have a transparent connection 
between the transformations of fields and those of the 
quantum states that they create. "Unbarred" operators 
generally involve the same constructions except that co- 
efficients or other operators involved must be hermitian 
conjugated. 

Operators that couple only to even parity interme- 
diate states in Eq. (11) are labeled with a subscript 
g (for gerade) and operators that couple only to odd 
parity states are labeled with a subscript u (for unger- 
ade). For half- integer spins, the relevant IRs of Of are: 
Gig, G2g 7 H gi Giu, Gi Ul H u . We will use these notations 
throughout this paper. 

Because of the parity transformation of Eq. (10), there 
are two independent ways to make baryon operators that 
couple to states of a given parity in a zero-momentum 
correlation function. Operators coupling to gerade states 
can be made either with even smearing parity together 
with positive p-parity or with odd smearing parity to- 
gether with negative p-parity. Similarly, there are two 
disjoint sets of operators that couple to ungerade states: 
ones with odd smearing parity together with positive p- 
parity or ones with even smearing parity together with 
negative p-parity. These sets are not connected by ro- 
tations because neither the smearing parity nor the p- 
parity can be changed by a rotation. However, they are 
connected by p-spin raising or lowering operations and 
in our construction each operator that couples to a ger- 
ade state is connected in this way with an operator that 
couples to an ungerade state. This is the second reason 
that the p, s labeling is useful. The p, s labeling is used 
sparingly in this paper but it is central to the method 
used in Appendix B to construct combinations of Dirac 
indices that transform irreducibly. 

Each baryon operator carries a row label, A, whose 
meaning depends upon the bases used for IRs. The row 
label distinguishes between the d\ members of IR A. If 
a representation contains more than one occurrence of 
IR A, we say that there are multiple embeddings of that 



IR. A superscript, k, is used to distinguish between the 
different embeddings. Therefore, a generic baryon oper- 

A k 

ator is denoted as B x ' (x, t), or in "unbarred" form as 
B^' k (x 7 t), where the operator belongs to the k th em- 
bedding of IR A and row A of the octahedral group. 
Operators for different baryons are indicated by the use 

A k 

of appropriate symbols, such as N \' (x, t) (for isospin 

A k 

1/2 operators), A A ' (x, t) (for isospin 3/2 operators), 

A k 

S A ' (x, t), and so on. 

The correlations of operators belonging to different IRs 
or to different rows of the same IR vanish: 

^(0|TB A A - fe (x,^'' fe '(0)|0) =C$(t)8 A A>6xx,. (15) 

X 

The correlations of different embeddings of the same IR 
and row are generally nonzero, providing sets of operators 
suitable for constructing a correlation matrix cj£} (t) . 

III. QUASI-LOCAL BARYONIC OPERATORS 

Since the quark fields are Grassmann- valued and taken 
at a common location x, and the color indices are con- 
tracted with the antisymmetric Levi-Civita tensor, our 
three-quark, quasi-local baryon operators must be sym- 
metric with respect to simultaneous exchange of flavor 
and Dirac indices. An operator that is symmetric in fla- 
vor labels (A, O) must be symmetric also in Dirac indices, 
and an operator that is mixed-antisymmetric in flavor la- 
bels (N) must be mixed antisymmetric in Dirac indices, 
assuming that masses of the up and down quarks are 
equal. An operator that is mixed-antisymmetric in fla- 
vor labels and that has nonzero strangeness (A) can have 
mixed-antisymmetric or totally antisymmetric Dirac in- 
dices, and an operator that is mixed-symmetric in flavor 
labels and that has nonzero strangeness (£, S) can have 
mixed-symmetric or totally symmetric Dirac indices. All 
possible symmetries of the Dirac indices are encountered 
in the consideration of the different baryons. In this sec- 
tion, we discuss the different baryons in turn and develop 
tables of operators classified according to IRs of . 

A. Quasi-local Nucleon Operators 

Consider operators made from quasi-local quark fields 
for isospin quantum numbers / = 1/2, I z = 1/2. These 
operators correspond to the N* family of baryons and 
they may be chosen to be 

where u is an up quark and d is a down quark. All 
(smeared) quark fields are defined at spacetime point 
(x,i). Equation (16) provides a proton operator in the 
notation of the Particle Data Group [37]. A neutron or 
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N operator can be obtained using the isospin lowering 
operation. 



e 



e 



MA 



MS 



FIG. 1: Young tableaux for three-quark 5(7(2)/ irreducible 
representations. The first tableau is antisymmetric in labels 
of particles I and 2 (denoted MA for mixed-antisymmetric) , 
while the second tableau is symmetric in the labels of particles 
1 and 2 (denoted MS for mixed-symmetric). The third tableau 
is fully symmetric (denoted S). 



These two operators correspond to the mixed- 
antisymmetric Young tableau for isospin in Fig. 1. 
Each iV$^jL 3 (x, t) operator of Eq. (16) is manifestly 
antisymmetric with respect to the flavor interchange 
u <-> d applied to the first two quark fields. This leads to 
the following restrictions on Dirac indices, 



N 



(MA) 

JilA»2At3 



N 



(MA) 



-(MA) -(MA) -(MA) 



= 0, 
= 0. 



(17) 
(18) 



Some general considerations are stated most simply 
using the Dirac indices. There are 4 3 =64 combinations 
of Dirac indices for operators formed from three quark 
fields. They may be classified by the four Young tableaux 
of Fig. 2, where each box is understood to take the values 
H = 1,2,3 or 4. 



S 
20 



e u e u e 

MS MA A 
20 20 4 



indices are written in accord with one Young tableau in 
Fig. 2 form a closed set under the group of rotations 
and parity transformations. These group representations 
have been reduced to Gi g / U and H g / U IRs of Dirac in- 
dices by working with the p, s labeling as discussed in 
Appendix B. 

Nucleon operators with MA isospin symmetry in 
Eq. (16) must have the MA symmetry of Dirac indices, 
corresponding to the third tableau of Fig. 2. Table VI 
gives explicit forms for the 20 quasi-local nucleon opera- 
tors classified into IRs G lg , H g , Gi u , and H u . Here and 
in the remainder of this paper we label the operators by 
the parity of intermediate states to which they couple 
in a zero-momentum correlation function, as in Eq. (11). 
Alternatively, one may regard the operators in the tables 
as having been translated to point x = 0, where they 
have definite parity as seen from Eq. (9). Dirac indices 
in the table come from Table XVI in Appendix B, but 
they have been simplified using the relation in Eq. (17). 
Because all the coefficients are real, "unbarred" operators 



TABLE VI: Quasi-local Nucleon operators. All operators 
have MA Dirac indices. 



^S,S Z "M1M2M3 


yA,fc -7jr= 
^S,S Z "M1M2M3 


tf? 1 ?' 1 N121 
F 2 1 

*i -i Nl22 


Vfr 1 -^(Nl23+Nl41+N 3 2l) 
*f'"l 7f(iVl24+iVl42+iV322) 


*r|' 7i 7|(iVl43+iV323+iV34l) 
71(^144+^324+^342) 


* ( r 1 i" 2 A343 
2 ' 2 

-=-G 1 u . 2 -t-t: 

1 1 JV344 

2 ' 2 


*f ^(A 7 13 4+A 7 323-A 7 34l) 

2, ■ 2 V d 

ft 2 VA — : 


*|'|' 3 7=(-Vl41-iV312+A 7 1 23) 

«? lu f 4- {-N322-N241+N124) 

2 ' 2 V3 


*3 9 3 V133 

h' 2 

*3 9 1 -=(^134+^143+^233) 

75(^144+^234+^243) 

*f 9 3 A244 
2 > 2 


*3"3 77l31 
2 ' 2 

*fl ^(iVl32+iVl41+iV 2 3l) 
*f"_l 73(^142+^232+^241) 

*f" 3 iV242 

2 ' 2 



FIG. 2: Young tableaux for three Dirac indices. 



Standard rules for counting the dimensions of the 
tableaux show that the totally symmetric tableau in- 
cludes 20 operators, the mixed-symmetric and mixed- 
antisymmetric tableaux each have 20 operators and the 
totally antisymmetric tableau has 4 operators, thus ac- 
counting for all 64 possibilities. Groupings of Dirac in- 
dices according to the symmetries of Fig. 2 are useful. 
The fact that a quasi-local baryonic operator must be 
symmetric with respect to simultaneous exchange of fla- 
vor labels and Dirac indices associates each baryon op- 
erator with one of the symmetries of Dirac indices found 
in the tableaux of Fig. 2. All the combinations of Dirac 
indices that correspond to each Young tableau are given 
explicitly in Appendix B. Three Dirac spinors whose spin 



are obtained by replacing TV by TV in the same linear com- 
binations . The left column of the table shows 10 gerade 
nucleon operators and the right one shows 10 ungerade 
operators. For a given parity there are three sets of G\ 
operators (three embeddings of G\) and one set of H op- 
erators. Each Gi IR contains two operators that trans- 
form amongst themselves under rotations of the group 
and each H IR contains four operators that transform 
amongst themselves. Operators in each IR are given spin 
projection labels, S z , which are also equivalent to "row" 
labels but more physically meaningful. In a given em- 
bedding the operator with the largest S z is designated 
"row 1", the next largest S z is designated "row 2", and 

so on. The notation 'Fg^ represents a general quasi- 
local baryonic operator with spin S and spin projection 
S z , transforming according to the k-ih embedding of IR 
A of the group . 
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Spin-raising and spin-lowering operators for a Dirac 
spinor are 



U\ ± i(7 2 

<7l ± if72 



(19) 



in the Dirac-Pauli representation. For a three-quark 
state, the spin raising or lowering operator is a sum of 
three terms, for example, S~ = X^j=i s f> where acts 
on the j-th quark. The same operations carry over to the 
"barred" field operators of Table VI. Different rows in 
the same embedding of an IR are related to one another 
by spin raising and lowering operations. For example, 
the transformation of the first G\ embedding of Table VI 
proceeds schematically as follows, 



-tttGi 



S-*v2,i/2 = S-N 121 =S-N+_ + 



= n 122 = *f; 2 9 ;! 1/2 ,(2o) 



Mixed-symmetric isospin operators with I = 1/2, I z = 
1/2 may also be defined by 

(23) 

However, for quasi-local operators they can be rewritten 
in terms of the MA isospin operators defined in Eq. (16) 
as follows, 



¥ (MS) =I (MA) -(MA) 



(24) 



showing that the quasi-local, MS isospin operators are 
not linearly independent of quasi-local MA operators. It 
is sufficient to consider only the MA operators of Eq. (16) 
in order to construct a complete, linearly independent set 
of isospin 1/2, quasi- local operators. 

Relations of the operators presented in Table VI to 
operators that are commonly used in the literature for a 
nucleon are given in Appendix C. 



where the notation AT^^f' 3 is used in the intermediate 
steps. Note that a spin-lowering operation on the sec- 
ond quark in Eq. (20) vanishes because it has spin down 
and spin-lowering of the first quark also vanishes because 
s^A^^ = N^t^ = by Eq. (17). Spin raising and 
lowering operations can be applied repeatedly and the 
following relation holds, 







w l/2,l/2 £jT *l/2,-l/2 



S+ 



S+ 

^3/2,3/2 W 3/2,l/2 



o. 



^3/2 -3/2 



H g/U ,k 
3/2,-1/2 

S 



0. (21) 



A gerade operator in a row of Table VI and the unger- 
ade one in the same row are related to each other by p- 
spin raising and lowering operations. For Dirac spinors, 
the operators are p ± = £\ =1 pf, pf = (l/2)(p] ± ipj), 
where j specifies the first, second, or the third quark. An 
example is 



P *V2,i/2 = P Ni2i=p n\ + _X, 

_L (n++- +N+-+ + N-++ 
^ ^ \ ly +-+ ^ ^ +-+ ^ ly +-+ 

-j-Gi„,l 
— *1/2,1/2j 



(22) 



where an appropriate normalization is included in the 
resultant operator. Note that p-spin raising and low- 
ering operations change the p-parity of one quark and 
thus change the product p\p 2 pz which is the p-parity of 
the operator. However, they preserve the s labels and 
leave the transformation properties under rotations un- 
changed. 



B. Quasi-local A and Q, Operators 

The isospin of a A baryon is 3/2 and there are four 
different operators corresponding to isospin projections 
I z = 3/2, 1/2, -1/2, and -3/2: 

A ++ 
A 



+ _ (abc l a h jc , -a lb -c , la -b -c \ 



+ ■■■ ■■ = 

- la lb 1c 
tabc « u , u u 2 u U3 ' 



-s-0 



~/\~ — r . Ja 3b Jc 



M1M2M3 



(25) 

where all fields are defined at spacetime point (x,i). Be- 
cause of the totally symmetric flavors, the A baryon oper- 
ators must have totally symmetric combinations of Dirac 
indices. According to Table XV in Appendix B there are 
20 combinations of totally symmetric Dirac indices. In a 
color-singlet three-quark operator, the quark fields may 
be commuted with one another with no change of sign. 
This allows Dirac indices to be rearranged to a standard 
order in which they do not decrease from left to right, 
producing the 20 irreducible operators that are given in 
Table VII. For each parity, two embeddings of the H IR 
occur, while there is one embedding of the G\ IR. Ta- 
ble VII holds for any I z value. Spin-raising and lowering 
operations as in Eq. (21) and p-spin raising and lowering 
operations as in Eq. (22) can be applied to relate opera- 
tors in different rows or operators in different columns of 
Table VII. 

The Q baryons are composed of three strange quarks, 



„ , ^ 3 ( x ' *) = 8°, (x, t)s b ^ (x, t)s% (x, t). (26) 



The quark flavors are clearly totally symmetric so only 
the totally symmetric Dirac indices are allowed. There- 
fore, Table VII can be used for an £1 baryon. In summary 
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TABLE VII: Quasi-local A operators. All operators have 
Dirac indices in S. 



TABLE VIII: Quasi-local A baryon operators. 



*S,S Z M1M2M3 


^S,S Z A M1M2M3 


^iY 1 A 13 4-A 23 3 
A144 — A234 


^I'l' 1 All4-Al23 

tff A124 - A223 

2 ' 2 


— tr~? 1 

* J fi 1 \/3 Sua 

*f 9 '\ V3A 122 
— Jr' ,i 2 — 

^3 _3 A222 


^tV \/3Aii3 

2 ' 2 

* f "1 1 A114 + 2A123 

2 ' 2 

^f"' 1 ! 2Al24 + A 22 3 

2 ' 2 

*f U,1 3 \/3A224 

2 ' 2 


¥3 9 3" \/3 A 133 
¥i 2 2 

* 3 9 ; 2 2A134 + A233 

*sfli Al44 + 2A 23 4 

*f 9 ' 2 3 V3A 2 44 

2 ' 2 


W 3 3 A333 

2 ' 2 

*f "1 2 Vs A334 

2 ' 2 

*?"' 2 i ^^344 
2 ' 2 

W 3 3 A444 

2 ' 2 



the A symbol in Table VII may be replaced with any of 
}■ 



{A ++ ,A + ,A°,A ,n 



W S,S Z M1M2M3 


V S,S Z il MlM2M3 


^V 1 !' 1 A 121 

2 ' 2 

W 1 1 iV122 


^T 1 ?' 1 A123 + A141 + A321 

2 ' 2 

\T/i* u \ A ln , 1 A1.ni A ...... 

*1 1 iV124 1 ^V142 T i\-322 

2 ' 2 


^f 1 ?' 2 A143 + A 323 + A341 

* 1 1 9 'i 2 A144 + A 324 + A342 


#1 1 A343 

2 ' 2 

*i 1 A344 
2 ' 2 


§ 

*l l9 ' A134 + A 3 23 ~ A 3 41 
g'2 3 

* 1 19 1 A144 + A 4 23 - A 3 42 


* i^r 1 -A141 - A312 + A123 
G 2 1 

\tl lu 'l — A322 — A 2 4i + A124 

2 ' 2 


* < 5 Jl r 4 a/|(Ai34+A 3 41+A413) 
2 ' 2 V 

*? l9 'l 4 \/f(A234+A 3 42+A 4 23) 


*f'|' 4 (Al23+A231+A 3 l 2 ) 
*? 1U 'l 4 -\/f(Ai24+A241+A 4 12) 

5.-5 V 3 


* 3 9 3 Al33 

— Jr' 2 — — — 

*3 9 1 A134 + Al43 + A233 

— Jr' 2 — — — 
3 1 A144 + A234 + A243 

— Jr' 2 — 

^3 3 A244 

2 f 2 


W3 3 A131 

2 ' 2 

* 3 "1 A132 + A 141 + A231 
H 2 

$ 3 1 Al42 + A 2 32 + A24I 

h' 2 

3 " 3 A242 

2 ' 2 



C. Quasi- local A Baryon Operators 



Irreducible basis operators for A c and A& baryons are 
exactly the same except that the third quark is replaced 
by a charm or bottom quark. 



The A baryons have isospin zero and strangeness — 1. 
Appropriate quasi-local A baryon operators have the 
form, 



1 



€abc \/2 ^ 2 ^ U ^ 2 ^ ^ 



(27) 



where spacetime arguments (x, t) are omitted from each 
quark field. The A baryon operator has a pair of up and 
down quarks in the isospin zero state, which is the same 
as the mixed-antisymmetric nucleon operator. Because 
the operator in Eq. (27) satisfies the relation 



^HlH2H3 ^A*2MlM3 



(28) 



it is antisymmetric with respect to exchange of [i\ and 
fj,2 indices. Allowed symmetries of Dirac indices for the 
quasi-local A baryon operator are mixed-antisymmetric 
and totally antisymmetric. The difference from the quasi- 
local nucleon operator is that the A baryon operator 
is allowed to have totally antisymmetric Dirac indices, 
because the strange quark removes the restriction of 
Eq. (18). Table VIII gives all quasi-local A baryon op- 
erators. Twelve positive-parity operators are given in 
the left half of the table and twelve negative-parity opera- 
tors are given in the right half. Only four combinations of 
Dirac indices are totally antisymmetric under exchange, 
and they belong to G\ IRs. Together with the three 
embeddings of G\ that come from mixed-antisymmetric 
combinations of Dirac indices, this provides a total of four 
embeddings of G\ in each parity plus one embedding of 
the H IR for quasi-local A baryon operators. 



D. Quasi- local S and S Operators 

A S baryon has two light quarks forming an isospin 
triplet combination and a strange quark. Suitable S op- 
erators arc defined such that the first two Dirac indices 
refer to the light quarks, 



M3 ' 



Such operators satisfy the relation, 



^A»1^2M3 ^(U2/il(Ii3 — 0, 



(29) 



(30) 



showing that the Dirac indices must be totally symmetric 
or mixed-symmetric. 

A S baryon has two strange quarks and one light quark 
forming an isospin doublet, 



—0 



6 abc s p, 1 s H2^fJ. 3 



(31) 



Again the operators are symmetric under the exchange 
of fii and fi2- Thus, the allowed combinations of Dirac 
indices are the same as for the quasi-local £ baryon op- 
erators. Table IX presents all operators with symmetric 
and mixed-symmetric Dirac indices. Note that there are 
20 operators for totally symmetric Dirac indices (as in 
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TABLE IX: Quasi-local E operators. 



A k 

^c'c S U1 „ oliq notation 


\ k 

'&<?'<? E U1U , u „ notation 


^l' 1 ?' "^[—^332 — 2S233 + S341 + E413 + £134] 

h' 2 1 

^ l 1S _'l -^[2Si44 + S441 — £234 — £342 — E423] 


* ' ^ [E114 + 2E141 — E123 — E312 — E231] 
a 2 1 

^ ~sr[ — ^223 — 2E232 + S124 + S24I + S412] 


tf? 1 !' 2 J||Sii a -Ei 2 i] 

2 ' 2 V 
T? 19 '! 2 \/|[-E 2 21 +S122] 

5. 2 V 


1 1 1 ' [E114 + 2E312 — E123 — E141 — E231] 

2 ' 2 a 

1 lu 'l \— E223 — 2E241 + E124 + E232 + E412] 

2 ' 2 3 


* 1 1 1 ' -^r [S332 + 2S134 — S341 — S233 — S413] 

2 ' 2 

*l l9 'l S441 — 2S423 + E342 + E144 + E234] 


^l 1 ?' \/lp334 — E343] 
2 ' 2 V d 

$1 y |[— E443 + E344] 


&' 1 

^1 V E332 ^ 2E413 + E233 + El34 + E34l] 

2 ' 2 

yrT^lgi^ r yT 1 q"v yT 1 
Wl 1 -7t-[2j441 + Z2j234 — 2^144 — 2^423 — 2^342] 


^l' 1 ?' ^[Ell4 + 2E231 — S141 — E312 — E123] 
2 ' 2 3 

TttGi u ,4 \/2 r Q 1 yT 1 I 1 
* A, -A 3 [ _ 2-.223 — 2Zj412 + 2j232 + 2j241 + 2j124J 


—ft £ " = 

^33 Em 

2 ' 2 

l^ 9 ; 1 73P112 +2E121] 
$3" 1 -75 [2E122 + E221] 

ft 1 2 

W 3 _ 3 2^222 


■^[En3+2E 13 i] 
^|"| |[2Ei4i + E114 + 2E312 + 2E123 + 2E231] 

* 3 1 Tj [2E232 + S223 + 2E412 + 2El24 + 2E241] 

2 ' 2 A 

vf7^ u ' * 1 1 1 

W3 _3 [2^224 + Z1j242\ 


*|J ^[2El33 + E33l] 

* 3 "i 77 [2E233 + E332 + 2El34 + 2E341 + 2E413] 

—ti 2 2 1 — — — — — 
*3 9 ' 1 tt[2Ei44 + E441 + 2E234 + 2E342 + 2E423] 

2 ' 2 

*f 9 ' 2 3 (l/v / 3)[2E 24 4 + E442] 

2-L—2 


T H„,2 7= 
W3 3 E333 
2 ' 2 

*f;i 2 -^[E334 + 2E 3 43] 

^S^ll [2E344 + E443] 
■J-H u ,2 

W 3 3 E444 
2 ' 2 


V|[-E331+El33] 

^ _ V _ _ _ 

' - S332 — 2E341 + E134 + E233 + E413] 
ft 2 3 

*3 9 ' 1 2E342 — E441 + E234 + E144 + E423] 

2 ' 2 

^ 3 a _3 ^/f"[— E442 + E244] 


*fr§ 3 ^/|[En3 - E131] 

\E , 3"i -r^[Ell4 + 2Ei23 — S312 — E141 — E231] 
H 2 

3 1 [2El24 + E223 — S412 — E232 — E241] 

2 ' 2 3 _ 

*|"J.| 1 [E224 — S242] 



the quasi-local A operators) and 20 operators for mixed- 
symmetric Dirac indices, giving a total of 40 operators 
for S or S baryons. Four G\ embeddings and three H 
embeddings occur in each parity. In Table IX the symbol 
£ may be replaced by any of {S , £ , £ , S , S }. 



IV. NONLOCAL BARYONIC OPERATORS 

In this section we discuss how to construct baryon op- 
erators that create states whose wave functions have an- 
gular or radial excitation. Orbital angular momentum or 
radial excitation is expected to be of particular interest 
for operators that couple to excited baryons. 

In Section III, all possible symmetries of the Dirac in- 
dices of three quarks were encountered. When nonlocal 
operators are constructed, we can build upon the quasi- 
local operators already found by adding a nontrivial spa- 
tial structure. This basically amounts to allowing differ- 
ent smearings of the quark fields. 

Nonlocal operators are constructed by displacing at 
least one quark from the others. The set of displacements 
is first arranged to belong to the basis of IRs of the oc- 



tahedral group. Then there arises the issue of combining 
the IRs of spatial distributions of displacements with IRs 
of the Dirac indices that have been developed for quasi- 
local operators. With respect to the octahedral group, 
the spatial and spin IRs transform as direct products. 
Using Clcbsch-Gordan coefficients, we form linear com- 
binations of the direct products so as to obtain nonlocal 
operators that transform as overall IRs of the group. 

A. Displaced quark fields and IRs of O 

Relative displacement of quarks requires insertion of 
a path-dependent gauge link in order to maintain gauge 
invariance. The simplest such displaced three-quark op- 
erator would be of the form, 

where the time argument is omitted from quark fields, 
and i is one of the six spatial directions {±x, ±y,±z}. 
Each quark field is smeared but the third quark has an 
extra displacement by one site from the other two in 
Eq. (32). 
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Spatial displacements of Eq. (32) with i <E 
{±x, ±y, ±z} transform amongst themselves under the 
rotations of the octahedral group O assuming that gauge 
links are cubically invariant (approximately true for av- 
erages over large sets of gauge- field configurations). The 
six-dimensional representation of O that is formed by the 
six displacements can be reduced to the IRs A\, T\ and 
E. In order to combine displacements so that they trans- 
form in the same way as the basis vectors of A\ , E, T\ 
IRs given in Table II, we first define the following even 
and odd combinations of forward and backward displace- 
ments: 



q ,Ji/ija _ rW/l/2/3 1 7"( 



fl Si /3 

fi h /3 



)hhh 



rW/1/2/3 T(-i)hhh 



(33) 



with i — x,y, z. The difference of forward and back- 
ward displacements, Di, has negative smearing-parity 
and involves a lattice first-derivative, while the sum of 
forward and backward displacements, Si, has positive 
parity. Note that the lattice first-derivative is an anti- 
hcrmitian operator. The second step is to form IR oper- 
ators using the Si and L>i combinations as follows: 



i 7-/1/2/3 / 



(x 



^ P1M2M3^ X 

T$u S ul ( X 

1 £il^2M3^ 



1 1 P2M3 y j 



Tfbi lh J 3 u (x 



1 

1 



E 



fi h /3 



(x), 



=1,2,3 



2\/3 \ 

A 7/1/2/3 

-o 2 o .„ 



S*i6 



/1/2/3 



(34) 

(x) 
(35) 



(x)), 

^(^™W-^^ 3 W) ,(36) 



2 (/>,*' M 



2 



A*lM2M3 

/l /2 /3 
H1H2H3 



(X) + iDyb^^i*) 



(x), 



(37) 



* A T/1/2/3 / \ 
7=E> z o,, ,, ,, x 

z H1H2H3 v / 



~2 



fx) 



(38) 



(^ 5 ™ 3 ( x ) 



• rS T/1/2/3 / \ 



A _ T/1/2/3 / \ 



(39) 



These definitions produce spatial distributions 
Ai,E*,T* that transform in the same way as the 
lattice harmonics of Table II. Superscripts on i? and Ti 
operators refer to the rows of the IRs. For the T\ ' ' 
combinations of displacements, we will generally denote 
operators by using the spherical notation as 
defined by Eqs. (37-39). 

We refer to these simplest nonlocal operators, involving 
linear combinations of operators with the third quark 
field displaced by one lattice site, as one-link operators. 
Let us denote the general form of a one-link operator as 



7-JAT/1/2/3 

A %ll^2/J3 



x)^ £abc ^{i(x)^ 



(x)[2?^ 3 (x)] c , (40) 



where A specifies the type of spatial IR (A\, T\ or E) 
and A specifies the row of the IR. In order to combine 
the spatial IRs of the displacement operators with the 
IRs of Dirac indices, we need the direct product rules. 



Direct products and Clebsch-Gordan 
coefficients 



Nonlocal operators involve direct products of two dif- 
ferent IRs of the octahedral group, one associated with 
the combinations of displacement operators and the other 
associated with the Dirac indices. Linear combinations 
of such direct products can be formed so that they trans- 
form irreducibly amongst themselves by using Clebsch- 
Gordan coefficients for the octahedral group. These have 
been published by Altmann and Herzig [21]. 

Clebsch-Gordan coefficients depend upon the basis of 
IR operators but different choices of the bases are re- 
lated to one another by unitary transformations. Be- 
cause our basis operators differ from those published by 
Altmann and Herzig, we have performed the required 
unitary transformations and obtained suitable Clebsch- 
Gordan coefficients for all possible direct products of 
IRs of the double octahedral group. A complete set of 
Clebsch-Gordan coefficients is given in Appendix D of 
this preprint and in Ref. [38] . The relative phases of op- 
erators from different rows within an IR should be fixed 
in lattice calculations in order to allow averaging over 
rows when that is appropriate, as it is in mass calcula- 
tions. However, different ways of forming a given IR as 
direct products need not have the same overall phases. 
We have used this freedom to eliminate phases within 
each table of Clebsch-Gordan coefficients such that all of 
our coefficients are real. 

A one-link operator that transforms as overall IR A 
and row A of D is written as a linear combination of 
displacement operators acting on IRs of Dirac indices as 
follows, 



= E c 



A1.A2 



A Ai A 2 
A Ai A2 



Ai A 2 



(41) 



where the corresponding quasi-local baryon operator is 

written as \I/ A ' instead of ^s's z an d the relation of A and 
S, S z is obvious from Table IV. For one-link operators, we 
need direct products of the IRs of displacements (Ai = 
Ai, E and T\) with the IRs of Dirac indices of quasi- local 
baryon operators (A 2 = G\ and H) . The following rules 
of group multiplication show which overall IRs A can be 
produced, 

A l ®G 1 = G u 
Ai ® H = H, 

E®d = H, 

E® H = G 1 (SG 2 ®H, 
TiOGi = G 1 ®H, 

Ti <8> H = G 1 ®G 2 ®H®H. (42) 
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C. One-link operators 

Baryon operators with one-link displacements can be 
categorized into two sets, one with antisymmetric and 
the other with symmetric Dirac indices of the first two 
quarks. The antisymmetric category includes the nucleon 
with MA isospin and the A baryon operators. The sym- 
metric category includes the nucleon with MS isospin, 
and the A, Q, £ and 5 baryon operators. These sym- 
metries determine the spinorial structures of the one- link 
operators. 

One-link operators for the nucleon with MA isospin 
and for the A baryon are taken to be of the form, 



_(3)A-T7(MA) _ e abc 
A JV MiA*2M3 ^2 

^abc 

7f 



u a d b - T u b 



V (3)A A 



m a*2 



Tl a —b 



where the superscript 3 of T>^ A denotes that the dis- 
placement operator T> A S {Ai, E x , T x } acts on the third 
quark, and 7V( MA ) denotes the isospin symmetry. This 
choice of one-link operators preserves the antisymmetry 
under /i! <-> /i 2 , and therefore requires Dirac indices to 
be MA (20 combinations) or A (4 combinations). Tak- 
ing into account the six possible T> A combinations of dis- 
placements, the total number of operators of the form of 
Eq. (43) or Eq. (44) is (20 + 4) x 6 = 144. 

One-link operators for the nucleon with MS isospin, or 
for the A, Q, £ and S baryons have the following forms: 



A(3)-^(MS) 

^\ iV Ml/J2M3 



€abc 



2« 









(45) 


D A(3) 5 ++ 

A ^1/^2^3 


= e abc U fj, 1 u fj, 2 




(46) 


U \ S VlP2M3 


—a —b 

— tabc S^S^ 




(47) 




= tabc u ^i l u ^2 




(48) 


T) A (3)^° 

A "^1/^2^3 


= tabc s fi 1 s i_i 2 


K% 3 ) c - 


(49) 



These operators are symmetric under fi± <-> /i2, so the 
allowed combinations of Dirac indices are totally sym- 
metric (20 combinations) or mixed-symmetric (20 com- 
binations). There are (20 + 20) x 6 = 240 such operators 
for each baryon. 



1. Ai one-link operators 

The reduction is the simplest for the A\ combination 
of one-link operators because it is just a scalar "smear- 
ing" . We show it as a first example. The MA isospin nu- 
cleon operator of Eq. (43) and the A baryon operator of 
Eq. (44) have the same restriction on Dirac indices as in 
Eq. (28). Because the A\ combination of displacements 
is cubically symmetric, these operators have the same 



transformations under group rotations as the quasi-local 
A baryon operators in Eq. (27), except that the strange 
quark is replaced by (Aiu^ 3 ) c and (Ais^Y, respectively. 

Dirac indices for A^ N and A are obtained from 
Table VIII, the quasi-local A baryon operator table. 

For each operator in Eqs. (45)-(49), the displacement 
makes the third quark distinct but the operators are sym- 
metric under [i\ «-> fi2 as in Eq. (30). This means that 
these operators transform in the same manner as the 
quasi-local £ baryon operators and Table IX can be used 
for any of the operators in Eqs. (45)-(49). 

We note in passing that any cubically symmetric form 
of smearing can be developed by repeated application of 
the A\ operator. Thus, any such smearing that makes the 
spatial distribution of the third quark different from that 
of the first two can be substituted for the A\ combination 
of displacements of the third quark. All such operators 
have the same transformations and thus the same IRs of 
Dirac indices. 



2. Ti,E one-link operators 

In order to construct operators that have the T\ or 
E combinations of one- link displacements, we apply the 
Clebsch-Gordan formula of Eq. (41) using the coefficients 
for the double octahedral group from Appendix D. The 
resulting one-link operators are given in Table X and Ta- 
ble XI. These tables give all possible T\ and E one-link 
baryon operators. The parity labels refer to the inter- 
mediate states that the operator couples with in a zero- 
momentum correlation function. In Table X we employ 
the notation of D + '~'°, instead of T 1 ' ' . The displace- 
ments are understood to act on the third quark. These 
tables are general in the sense that they apply to any 
baryon, e.g., N, A, A, E, or H. 

A k 

The notation ^ x describes a quasi-local operator 
whose spin belongs to the fc-th embedding of IR A and 
row A. These operators are taken directly from the ta- 
bles for quasi-local baryon operators discussed in Sec- 
tion III in a similar fashion as for the A\ one-link oper- 
ators. One-link nucleon operators with MA isospin and 
one-link A baryon operators employ the spinorial struc- 

tures ^ x used for quasi-local A baryon operators given 
in Table VIII, together with Table X for 7\ and Table XI 
for E. One-link nucleon operators with MS isospin or 
one-link A, Q £ and S operators employ the same spino- 

A k 

rial structures \& A ' that are given for the quasi-local £ 
baryon operators in Table IX. 

The Ti one-link operators in Table X are strictly 
"barred" fields. The corresponding "unbarred" T\ op- 
erators use hermitian conjugated lattice first-derivatives. 
The factor i that has been included in the T± one-link op- 
erators provides the same hermiticity property as spher- 
ical harmonics, i.e., Y l \ n = (— l) m Yj ; _ TO . Note that be- 
cause the smearing parity of the T\ displacement is nega- 
tive, the overall p-parity is opposite to the overall parity 
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TABLE X: Ti one- link operators. Note that D + = fl, D° = Ti 2 , and D~ =ff. 



IR row C( A Ti A 2 ) J D(3)^i*^' fc 

\ A Ai A2/ 


- |J, J,> 


h 9/u i d^T'* 

2 '. 2 -G, fe 2 ' 2 


2 + 3\ 

2 ,I\ 
2' T 2/ 

3 _i\ 

2 ' 2 / 

3 3 \ 
2 ' 2/ 


V 2-2 » 2 ' 2 

V - r > *1 3 T ^[f t U *3 1 A/ 15 *3 _1 

V 2'2 V1J 2'2 V 2'2 

V° 2 i 2 V° 2 > 2 


3 ,3\ 
2 |T 2/ 

3 ,I\ 
2' T 2/ 

3 _I\ 

2 ' 2 / 

3 3\ 
2 ' 2 / 


2 -^o^-^fn 3 '" - ^°^fr/ s ' fc - v^¥^ + *f^i' fc i§.+*> 

d V 10 §,i + V 5 f,-i + v^o W f,-| la- 2/ 


V 2 ' 2 2'2 

2 1 EPtf?»';' k -Jib- ^J" k 


2 |T 2/ 
2 ' 2 / 


2 ^ *|i V3 W +72^ 


- +-) 

2' T 2/ 

2 ' 2 / 


r> 1 1 n-W H " /3 ' fe 1 n°vi7 H " /9 '' i: -l 1 n+vfr^"/ 9 '* 1 I 5 /s~|s 3 \ 

G 29/U 1 -75O *3_1 -7f£>* 2i _3 *3 3 75I2.+2)- V 6 l3'-2> 
2 ^ U W §_3 *3 3 ^ V 6 l2' + 2/ + | 2 > 2/ 



provide two members of the rank-two spherical harmon- 
ics. The remaining three members belong to the T2 IR 
and they cannot be constructed unless there are least 
two displacements in perpendicular directions, as will be 
discussed in the next section. 

For baryon fields with projection to zero total momen- 
tum, the following linear dependence holds, 

x 

^^(x)(/>% 2 (x))^ 3 (x) 

X 

- ^(^% i (x))g Al2 (x)^ 3 (x) + 0(a 2 ). (50) 

X 

This relation derives from the fact that after projection 
to zero total momentum, a total derivative of a baryon 
field vanishes and a total derivative is equivalent to or- 
der a 2 to a sum of lattice derivatives applied to each 
quark field. Some of the T± one-link nucleon operators 
are not linearly independent because of this. A nucleon 
operator with MS isospin having MS Dirac indices is 
equivalent (within a total derivative) to a nucleon op- 
erator with MA isospin having MA Dirac indices, for 



for T\ one- link operators. 

The last column of Table X shows the total angular 
momentum and its projection onto the z-axis. The first 
set of H operators involves the same Clcbsch-Gordan co- 
efficients as would apply to the formation of continuum 
states: (L = 1) ® (S = 1/2) — > (J = 3/2). This is a 
result of using the combinations of displacements that 
transform in the same way as the basis vectors of Ta- 
ble II. The second set of H operators has the continuum 
Clebsch-Gordan coefficients for (L = 1) ® (S = 3/2) — ► 
(J = 3/2). The third set of H operators corresponds to 
(L = 1) ® (S = 3/2) — > (J = 5/2), but J z values are 
mixed in row 1 and row 4. Similarly, the first set of G\ 
operators in Table X is (L = 1)®(S = 1/2) -> (J = 1/2), 
and the second one is (L = 1) ® (S = 3/2) — > (J = 1/2), 
both having the continuum Clebsch-Gordan coefficients. 
The G2 operators correspond to (L = 1) ® (S = 3/2) — > 
(J = 5/2), but J z values are mixed. 

Direct products involving the E spatial IR of dis- 
placements and spinorial IRs are given in Table XI. 
No operators involve continuum Clebsch-Gordan coeffi- 
cients in Table XI because the E IR has mixed J z , i.e., 
E 2 ~ F2,2 + 5"2-2- The E combinations of displacements 
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TABLE XI: E one-link operators. All operators have mixed 



TT> 

IK 


row 




H g / U 


1 


CI h 




2 


2 ' 2 
2 ' 2 




3 


-E^^'r" 




4 


2 ' 2 


Hg/u 


1 


- 1 — ff / k - o — ff / fc 

-s 1 *^ 1 " - £ 2 *3 9/ T 




2 


2 ' 2 2 ' 2 

E 1 * 3 fl ( u - E 2 V 3 S/ 3 




3 


2 ' 2 2 ' 2 




4 


2 ' 2 2 ' 2 

— E ^ 3 9/ 3 — £/ ^31 
9.'- 9. 9,19. 


Glg/u 


1 


E^lf k + E^f'_f 




2 


-E^f^ k -E^ H /t k 

9. • 9 9.' 9. 


G'lg/u 


1 






2 


9' 9 9' 9 



TABLE XII: Allowed combinations of Dirac indices for differ- 
ent one- link (Ai,Ti,E) baryons. The displacement is always 
taken on the third quark for simplicity. The third quark of 
the A and S baryons is chosen to be strange quark, and the 
third quark of the H baryon is the light quark. The numbers 
of operators for A\, Ti, or E combinations of displacements 
are listed in the fifth column and the numbers of operators 
for each overall IR of O d are shown in the last three columns, 
counting both parities, all embeddings and all rows. Linear 
dependencies resulting from a projection to zero total mo- 
mentum are not taken into account in this table. 



one-link baryon 


spin sym. 


Table 


disp. 


# 


d 


G 2 


H 


iV< MA ),A 


MA, A 


VIII 


Ai 


24 


16 





8 






X, VIII 


Ti 


72 


20 


4 


48 






XI, VIII 


E 


48 


4 


4 


40 


jV( MS ),A,n,E,H 


MS, S 


IX 


A-l 


40 


16 





24 






X,IX 


Ti 


120 


28 


12 


80 






X,IX 


E 


80 


12 


12 


56 



follows, 



the T\ one-link construction. It is easy to show that 
Di 3) N { ™2v3 + 1(^2 " Ms) + Ma) can be writ- 

ten as a linear combination of J s by applying 

Eq. (50). This identity reduces the number of Ti one-link 
nucleon operators by 20 x 3 = 60, where the number 20 
comes from the number of MS Dirac indices (or MA Dirac 
indices) . The number of distinct one-link nucleon opera- 
tors (both MA and MS isospin) after projection to zero 
total momentum is 64 for A\, 132 for Ti, and 128 for E. 
The total number is 324. 

Operators that are totally symmetric with respect to 
flavor exchanges, such as the A baryon, have a similar 
restriction. Such operators vanish when a first-derivative 
acts on one quark in a totally symmetric combination of 
Dirac indices. There are sixty A baryon operators with 
the Ti one- link displacements that vanish after projection 
to zero total momentum. 

The correspondence between the type of baryon and 
the symmetry of Dirac indices for the two categories of 
one-link baryon operators is summarized in Table XII. 
The numbers of possible operators are shown for con- 
structions using A\, Ti, or E spatial IRs to obtain G\, 
G2, or H overall IRs. 



T) (3)A 2T) (3)Ait/i/2/3 

U \2 U \l C W2M3 



(x) 



^(x)g^(x) 



(51) 



where the third quark is displaced covariantly by two 
displacement operators V^,V^ e {A 1 ,E X ,T^}. The 
first displacement acts on the third quark and defines a 



modified quark field, Q^f 3 ^ 1 (x) = 



X> Al 77^ 3 



(x) 



at the 



same position x. Then the second displacement further 
displaces the field and so defines a second modified field 

~c/ 3 ,Ai,A 2 . . 



at the same position, 

Figure 3 shows schematic illustrations of three distinct 
displacement configurations for a two-link baryon oper- 
ator. The first figure shows the "bent-link" operator, 
where a line denotes the gauge link and the arrow speci- 
fies the point at which the displaced quark's color index 
forms a color singlet with the other quarks. The sec- 
ond figure shows the possibility that the third quark is 
translated back to its original position by the second dis- 
placement, which is equivalent to a quasi- local operator 
because U l {x)u}{x) = 1. The third fi gure shows the pos- 
sibility of two displacements along the same direction, 
which gives a straight path differing from a one-link dis- 
placement only by its length. Inclusion of the bent-links 



D. Two-link operators 

One- link operators make it possible to realize A\, Ti 
and E types of spatial smearing, but not the T 2 or A 2 
types. The latter two types appear in the two-link op- 
erator constructions. We define a two-link operator as 



FIG. 3: Illustration of two-link displacements, 
can enrich the angular distribution and recover parts of 
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the continuum spherical harmonics that cannot be ob- 
tained from one-link displacements. 

First, we classify the spatial degrees of freedom into 
a single IR of O by forming linear combinations of the 
elemental operators of Eq. (51). The overall spatial IR 
A and row A are determined by the direct product of the 
two spatial displacements T>^ and with appropriate 
Clebsch-Gordan coefficients, 



E c 

Ai ,A 2 



A Aj 
A Ai 



A2 Ai 



(52) 



A particular example is instructive. Suppose one chooses 
V^l to belong to the 7\ IR, and X>^ 2 to belong to the 
E IR and desires the overall spatial IR to be T\. Then 
Eq. (52) is used with Clebsch-Gordan coefficients from the 
E ® Ti table in Appendix D, which gives 

^-E 2 x if, 
2 n 





2 


x fl 


fi ~ 


-E 1 


x fl 




2 E 


x fl 



y^-E 2 x fi. 

2 1 

In this way the two-link operator is determined so that 
its spatial part transforms according to a particular IR 
(in this case Ti). Once the overall spatial IR is obtained, 
Clebsch-Gordan coefficients for the direct products of the 
overall spatial IR and a selected spinorial IR are used to 
form an operator that overall transforms irreducibly ac- 
cording to Gi , H , or G2 • Because the spatial IR in the ex- 
ample above is Ti, which has been considered already in 
the construction of one-link operators, Table X provides 
the result. The only change is to use the T\ distribution 
of two-link displacements in place of the T\ distribution 
of one-link displacements. The use of Clebsch-Gordan 
coefficients of the cubic group has reduced the problem 
of finding IRs of two- link baryon operators to the already 
solved problem of finding one-link baryon operators. 

However, new possibilities exist with the two-link dis- 
placements. One can form the Ai and T 2 spatial IRs 
that did not appear in the one-link construction. This 
construction is straightforward but is omitted from this 
paper except to note that two-link A 2 and T 2 lattice har- 
monics correspond to the spherical harmonics shown in 
Table II. 

Proceeding in this fashion, one may construct multi- 
link baryon operators 



V 



(3)A„ 



(3)A 2 (3)^/2/3 ( } 
A 2 Ai Ml/ J 2M3 1 > / 



ea6c^ 1 (x)^(x) 



,(53) 



that involve n-site displacements in space allowing a 
quark to be displaced over a finer angular distribution so 
as to yield higher rank spherical harmonics. The reduc- 
tion procedure is essentially the same as for the two-link 
case, except that multiple direct products of spatial IRs 
are used. 



E. One-link displacements applied to two different 
quarks 

Consider an operator with one-link displacements ap- 
plied to two different quarks in the following way, 

t .(1±2)A 2t) (3)Ait/i/2/3 _ 

U X 2 U \l °MlM2P3 — 

€abc 

7f . 



where T> 



(1±2)A 
A 



indicates that the first two quark fields are 



symmetric or antisymmetric with respect to exchange of 
their spatial dependencies, 



r) (i±2)A _ 1 

~7i 



)A 



±Pf A ) 



(55) 



We refer to £)^ 1+2 ) A as space- symmetric and to 2 ^ A as 
space- antisymmetric combinations of displacements. The 
symmetry of the spatial displacements must be taken into 
account in the overall antisymmetry of operators in order 
to identify the symmetry of Dirac indices that produces 
nonvanishing operators. 

For the case of MA isospin nucleon operators with one- 
link displacements applied to two quarks, we obtain 

T) (1±2)A 2 „(3)A 1 ^t(MA) 

A2 Ml M2 A 2 A 1 ! 



.4 K^A -v a a<) K^y 

»4 -^>a 2 2 <) (<^ 3 ) C ,(56) 

and the following relation between spatial symmetry of 
displacements and the symmetry of Dirac indices holds 



± e 



T ,(l+2)A ! „(3)A 1 Tj(MA) 
^A 2 ^Ai J >iA»2M3 



A, MA Dirac indices 



V 



(l-2)A 2T .(3)A 1 T f (MA) 



^r 1 ^™ - s > ms ° irac indices - ( 57 ) 



Because of symmetry, the operators of Eq. (56) with MA 
Dirac indices and those with MS Dirac indices are iden- 
tical. 

Group theoretically, rotations of operators with one- 
link displacements applied to two different quarks are the 
same as those of two-link operators. Therefore the reduc- 
tion to IRs is exactly the same as for the two-link case. 
First use the Clebsch-Gordan coefficients to obtain an IR 
for the product of two displacements, and second use the 
Clebsch-Gordan coefficients for the direct product of spa- 
tial and spinorial IRs to obtain operators corresponding 
to overall IRs. The only additional step is to determine 
the allowed symmetries of Dirac indices such that the op- 
erator is antisymmetric under simultaneous exchange of 
displacements, flavors, colors, and Dirac indices. 



V. SUMMARY 

The constructions given in this paper provide a variety 
of quasi-local and nonlocal three-quark operators for use 
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as zero-momentum baryon interpolating field operators 
in lattice QCD simulations. All operators are catego- 
rized into the double- valued IRs of the octahedral group 
O®, they have definite parities and they are gauge invari- 
ant. Operators correspond as closely as possible to the 
continuum | Jm) IRs and they should be useful for spec- 
troscopy and for applications that require baryons with 
a definite spin projection. 

Complete sets of quasi-local operators are presented 
in Section III for each baryon. These quasi-local con- 
structions provide templates for the Dirac indices that 
should be used to construct nonlocal operators. Nonlocal 
operators are developed in Section IV based on adding 
combinations of one-link displacements to one or more 
quarks. By use of the building blocks given in this paper, 
a variety of additional operators can be constructed by 
1) using the Clebsch-Gordan series to form overall IRs of 
the spatial distribution, and 2) combining the spatial IRs 
with IRs of Dirac indices to form operators corresponding 
to overall IRs. Identification of the correct symmetry of 
Dirac indices is straightforward when space-symmetric or 
space-antisymmetric combinations of displacements are 
used. 

Reference [4] has demonstrated numerically that our 
quasi-local and one-link operators are orthogonal in the 
sense of Eq. (15), i.e., a correlation function vanishes if 
sink and source operators belong to different IRs and 
rows. For calculations of baryon masses, one should se- 
lect source operators within a fixed IR and row from the 
various tables. Using operators from different embed- 
dings of the IR and row, matrices of correlation functions 
may be calculated and mass spectra extracted. Correla- 
tion matrices can be made hermitian by including a 74 

matrix for each quark in the source operator. Operators 

\ k 

from our tables have the form B ' 



„(A,A,fc)-5 



where B,, 



is an elemental baryon operator and a 



summation over repeated indices is understood. A her- 
mitian matrix of correlation functions can be calculated 
in following way, 

C$ (*) - E <°l TB x k ^ t)B A x ' k ' (0) |o) r 4 



(A,A,fc)*„(A,A,fe') 



C t 1 



I TB. 



(x,t)fl 



-<,^3(0)|0) 

(58) 



Exploratory calculations for baryon spectra along this 
line have been reported in Ref. [3]. For a given baryon, 
the dimension of the matrix of correlation functions de- 
pends on the choices that are made for spatial distribu- 
tions (quasi-local, one- link, two- link, etc.) and the overall 
IR. For nucleon operators with quasi-local and one-link 
displacements, 23 G\ g operators, 28 H g operators, and 7 
G2 g operators are available as shown in Table XIII. The 
numbers of operators in each IR and row can be extended 
without limit by using two-link and three-link operators 
and by using different choices of smearing. 



TABLE XIII: Available numbers of nucleon operators with 
quasi- local and with one- link displacements for Gi g row 1, 
G23 row 1, and H g row 1. The numbers of ungerade operators 

are exactly the same. The Ti iV (MS) operators are omitted 
because after a projection to zero total momentum they are 
equivalent to the Ti N^ MA ^ operators. 



Type 


Eq. 


Table 


Gig 


H g 


G 2g 


quasi-local 


16 


VI 


3 


1 





one- link A\ 


43 


VIII 


4 


1 





one-link E 


43 


XI, VIII 


1 


5 


1 


one-link Ti 


43 


X, VIII 


5 


6 


1 


one- link A\ 


45 


IX 


4 


3 





one- link E 


45 


XI, IX 


3 


7 


3 


one-link Ti 


45 


X, IX 


3 


5 


2 


total 






23 


28 


7 
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A quark field expressed in terms of the Dirac-Pauli repre- 
sentation may be re-expressed in the Degrand- Rossi con- 
vention, for example, by 



g(DP) 



(V (DR)t ) 



(DR) 



(A5) 



In order to display the spin and p-parity of fields in a 
transparent way, we employ spin s subscripts and p su- 
perscripts in place of the four Dirac components p, = 1, 
2, 3 and 4 of each quark field in the Dirac-Pauli repre- 
sentation as shown in Table V. 

This encoding of the Dirac indices is based on the 
SU(2) ® SU(2) representation of the Dirac matrices, 
where the first SU{2) is generated by 2x2 Pauli matrices 
for p-spin, 



Pi = 




P2 = 







,P3 = 





-1 



(A6) 



and the second SU (2) is generated by the Pauli matrices 
for ordinary spin, 



ci = 




,02 = 



-i 

1 



,03 




.(A7) 



The 2x2 identity matrices for p-spin and ordinary spin 
are 



Pi = 04 




(A8) 



In terms of these sets of 2x2 matrices, the 4x4 Dirac- 
Pauli matrices are expressed as direct products of p-spin 
matrices and s-spin matrices as follows, 







J 


= Pi ® 


04, 


APPENDIX A: DIRAC MATRICES 












74 


= A3 ® 


04, 


Various conventions for the Dirac matrices are 


useful. 


75 


= Pi ® 


04, 


Each is related by a unitary transformation to the 


Dirac- 


7574 


= -ip2 


® 04, 


Pauli representation as follows, 




757fe 


= ipzd 


5 0fc, 




(Al) 


Ik 


= P2 ® 


Ofc, 


where 




04fc 


= Pi <g> 


Ofc, 






CfcZ 




m.Pi ® <J 



(DP) 
7 4 




(A9) 



(A2) 



where k, I, and m take the values 1, 2 and 3. 



The unitary matrix that generates the Weyl convention 
is 



, . (DP) (DP) 
Tj(W) = 1 + 7 5 74 

V2 



(A3) 



and the unitary transformation that generates the 
DeGrand- Rossi convention is 



, , _,^(DP) , ^(DP) (DP) 
jy(DR) _ 1 12 + 7l 73 

y/2 



(A4) 



APPENDIX B: SYMMETRY OF THREE DIRAC 
FIELDS 

The Dirac indices categorized in each Young tableau 
in Fig. 2 should be further reduced into G\ or H IRs 
for the purpose of operator construction. (There is no 
G 2 IR with three Dirac spinors.) Decomposition of the 
Dirac index into p-parity and ordinary two-component 
spin (s-spin) simplifies this process. 
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The S (totally symmetric), MS (mixed-symmetric), 
and MA (mixed-antisymmetric) combinations of three s- 
spins are defined as follows. 



S : |+++); 



!++-> + !+-+> + !-++>. 



V3 



l+— ) + !-+-> + 1— +> 



V3 



MS 



MA 



V6 



(2 - !+-+>- !-++»; 



--L (2 |- + )-|- + -)-l + -)). 



(Bl) 

(B2) 
-»(B3) 



The four states in Eq. (Bl) are \J,m) = |||), |||), 
|| — |), and ||— |) respectively. The two states in 
Eq. (B2) are \\\) and || — |) while the two states in 
Eq. (B3) are also \ \\) and || — i). All these states are 
orthogonal to one other. Because S states in Eq. (Bl) 
span total spin 3/2, they are the bases of an H IR (no 
matter which p's are involved in making up the Dirac 
indices). The MS and MA states in Eq. (B2, B3) span 
total spin 1/2, so they are the bases of G\ IRs. 

Products of three p-spins are categorized in exactly the 
same way. The p-parity is given by the product pip2P3- 
Direct products of states of three p-spins and states of 
three s-spins are simple when they are expressed in the 
bases of S, MS, and MA. For instance, MA p <g)S s with sub- 
scripts denoting p-spin and s-spin describes eight states, 
four of which have positive p-parity, and four of which 
have negative p-parity. The four states of each p-parity 
span H IRs because IRs of O d are determined only by the 
s-spins. The direct product of MA P ® S s , with ||,— 5) 

and || 5 1) is written as 



\+—) P - h+-) p ^ |++-) fl + |+-+) fl + |-++) s 



V2 V3 
By evaluating the direct product one obtains 

^ iWv-) + 11=;) + fe> - - \A) - 14;)) 



where the notation 



is used. 
IM1M2M3) notation, it becomes 



In terms of the 



1 



— (|134) + 1 143) + |233) - |314) - |323) - |413>) , 
V6 

where the Dirac indices p^ = 1,2,3,4 are defined in 
Dirac-Pauli representation of Dirac 7 matrices. The 
translation of p to (p, s) is given in Table V. It is clear 
that the obtained Dirac indices are antisymmetric under 
exchange of first two labels but not totally antisymmet- 
ric. Thus, we denote MA p ® S s = MAnirac- The nucleon 
operator that follows from this example is labeled as H g , 
row 2 in Table VI. 



From such considerations one obtains Table XIV, 
which provides the relations of Dirac index symmetries 
(abbreviated as "Dirac sym" in the table) to IRs of 
Dirac indices, and direct products of p-spins and s-spins. 
Note that MA„ <g> VIA, and MS„ <g> MS, both have a mix- 



TABLE XIV: Relation among Dirac spin symmetry, IR, and 
the direct product of p-spins and s-spins. 



Dirac sym 


IR 


emb 




p®s 




SDirac 


Gi 


1 


MA P $ 


5 MA S MSp C 


9MS 3 




H 


1,2 




Sp® S s 




MSDirac 


Gi 


1,2 




S p ® MS S 








3 


MA P <i 


5 MA S e MSp C 






H 


1 




MSp ® S s 




MA D irac 


Gi 


1,2 




Sp ® MA S 








3 


MA P <i 


5 MS S MSp $ 


}MA, 




H 


1 




MAp (g) S s 






Gi 


1,2 


MAp $ 


5 MS S e MSp $ 


}MA, 



ture of Soirac and MSoirac- One can easily see that addi- 
tion of a state from MA p ® MA S , say Gi g , row 1, and a 
state from MS p £g)MS s of the same Gi g , row 1 yields a pure 
SDirac state. The subtraction of the states yields a pure 
MS D irac state. Similarly, MA p ® MS S and MS P ® MA S 



have a mixture of MAr 



and A 



Dij 



A pure MAoi: 



is obtained by addition of states from MA p ® MS S and 
MS P (8) MA S and a pure Aoirac is obtained by subtrac- 
tion of states from MA p <g> MS S and MS P (& MA S . The 
third column of Table XIV shows an embedding that has 
a connection to Table XV in a self-explanatory way. 

Explicit combinations of Dirac indices P1P2P3 are 
given in Tables XV and XVI. Table XV contains all 
S and MS combinations of three Dirac indices, assigning 
each to an IR ( G\ g i u or H g / U ), embedding, and row. 
Table XVI contains all MA and A combinations of three 
Dirac indices in a similar way. 



APPENDIX C: RELATIONS OF V M1M2M3 TO 
COMMONLY USED NUCLEON OPERATORS 

Various groups have performed lattice simulations us- 
ing the following two interpolating fields for a nucleon: 



1/2 
Xi 

1/2 
X2 



u, 



(u T Cd) 75 u, 



(CI) 
(C2) 



where spacetime arguments are omitted. Matrix C is a 
charge-conjugation operator, defined by C — 7472- Each 
of these four-component operators corresponds to a G\ 

IR and may be written in terms of $ s ' s . Positive and 

1/2 

negative p-parity parts of \\ are projected in Dirac- 
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TABLE XV: IRs of S, MS combinations of Dirac indices for 
three quark states in Dirac-Pauli representation. The first 
entry of MS table reads 2q{ 1 q{ 2 q^ i -qtq^ql 3 -q^q^ql 3 for 
the Gi 9 , embedding 1, row 1 local operator. 



A 


k 


A 


S Dirac indices /^i^i2M3 


Gig 


1 


1 


-2(332 + 323 + 233) + 341 + 431 + 314 + 413 + 134 + 143 






2 


2(144 + 414 + 441) - 234 - 342 - 423 - 243 - 324 - 432 


G lu 


1 


1 


2(114 + 141 + 411) - 123 - 213 - 132 - 231 - 321 - 312 






2 


-2(223 + 232 + 322) + 214 + 124 + 241 + 142 + 421 + 412 


H g 


1 


1 


111 






2 


112 + 121 + 211 






3 


122 + 212 + 221 






4 


222 




2 


1 


133 + 313 + 331 






2 


233 + 323 + 332 + 134 + 341 + 413 + 143 + 431 + 314 






3 


t A A 1 ATA 1 A A 1 1 OQ/1 1 Q/IO 1 /I T Q 1 H A 1 1 1 *30/l 

144 + 414 + 441 + zo4 + o4z + 4zo + Z4o + 4oz + OZ4 






4 


244 + 424 + 442 


H u 


1 


1 


113 + 131 + 311 






2 


411 + 141 + 114 + 312 + 123 + 231 + 321 + 213 + 132 






3 


322 + 232 + 223 + 412 + 124 + 241 + 421 + 214 + 142 






4 


224 + 242 + 422 




2 


1 


333 






2 


334 + 343 + 433 






3 


344 + 434 + 443 






4 


444 


A 


k 


A 


MS Dirac indices fiifi2f^3 


Gig 


1 


1 


2(112) - 121 - 211 






2 


-2(221) + 212 + 122 




2 


1 


2(332 + 314 + 134) - 341 - 323 - 143 - 431 - 413 - 233 






2 


-2(441 + 423 + 243) + 432 + 414 + 234 + 342 + 324 + 144 




3 


1 


-2(332 + 413 + 143) + 323 + 233 + 134 + 314 + 341 + 431 






2 


2(441 + 324 + 234) - 414 - 144 - 243 - 423 - 432 - 342 


Gi„ 


1 


1 


2(114 + 132 + 312) - 123 - 141 - 321 - 213 - 231 - 411 






2 


-2(223 + 241 + 421) + 214 + 232 + 412 + 124 + 142 + 322 




2 


1 


2(334) - 343 - 433 






2 


-2(443) + 434 + 344 




3 


1 


2(114 + 231 + 321) - 141 - 411 - 312 - 132 - 123 - 213 






2 


-2(223 + 142 + 412) + 232 + 322 + 421 + 241 + 214 + 124 


H g 


1 


1 


-2(331) + 313 + 133 






2 


-2(332 + 341 + 431) + 314 + 134 + 323 + 143 + 413 + 233 






3 


-2(342 + 432 + 441) + 324 + 144 + 414 + 234 + 423 + 243 






4 


-2(442) + 424 + 244 




1 


1 


2(113) - 131 - 311 






2 


2(114 + 123 + 213) - 132 - 312 - 141 - 321 - 231 - 411 






3 


2(124 + 214 + 223) - 142 - 322 - 232 - 412 - 241 - 421 






4 


2(224) - 242 - 422 



TABLE XVI: IRs of MA, A combinations of Dirac indices for 
three-quark states in Dirac-Pauli representation. The caption 
in Table XV describes how to read this table. 



A 


k 


A 


MA Dirac indices 




Gig 


1 


1 


121 - 211 








2 


122 - 212 






2 


1 


143 - 233 + 323 - 413 + 341 - 


431 






2 


144 - 234 + 324 - 414 + 342 - 


432 




3 


1 


-233 + 323 + 134 - 314 - 341 - 


-431 






2 


-414 + 144 - 243 + 423 + 432 - 


- 342 




1 


1 


123 - 213 + 141 - 231 + 321 - 


411 






2 


124 - 214 + 142 - 232 + 322 - 


412 




2 


1 


343 - 433 








2 


344 - 434 






3 


1 


-141 + 411 - 312 + 132 + 123 - 


- 213 






2 


232 - 322 + 421 - 241 - 214 + 


124 


H g 


1 


1 


133 — 313 








2 


1 ^4 — ^1 4 4- 1 A^ — "V? 1 * 4- 9^3 


413 






Q 

o 


144 ^94 4-9^4 414 4-94*3 
I'll — OjZ*± 1 ZO^± — *±1*± "T i^io — 


423 






4 


244 - 424 






I 




131 - 311 








2 


132 - 312 + 141 - 321 + 231 - 


411 






3 


142 - 322 + 232 - 412 + 241 - 


421 






4 


242 - 422 




A 


k 


A 


A Dirac indices ^1^2/^3 




Gig 


1 


1 


134 - 314 + 341 - 431 + 413 - 


143 






2 


234 - 324 + 342 - 432 + 423 - 


243 


Glu 


1 


1 


-123 + 213 - 231 + 321 - 312 + 132 






2 


-124 + 214 - 241 + 421 - 412 + 142 



parity, 

1 + 74 1/2 = ( Nl43-^233 \ 
2 2 V iVl44 - ^234 J 

1-74 1/2 = I N U1 - N 23 i \ 
2 X2 ' \N 142 -N 232 J 

(C4) 

1/2 1/2 

These results show how the components of \\ an d X2 
are related to operators defined in this paper. 



Pauli representation as follows, 
1 + 74 1/2 

—z— Xi = 



1 - 74 1/2 

2 Xl ~ 



N 12 i 


-4V341 


N 122 


"^342 


N 123 


-^343 


N 124 


-^344 



(C3) 
f-1/2. Sim- 
ilarly x 2 ^ can be projected to operators of definite p- 



The upper component corresponds to S z 
1/2 



APPENDIX D: CLEBSCH-GORDAN 
COEFFICIENTS FOR THE DOUBLE 
OCTAHEDRAL GROUP 



The Clebsch-Gordan formula shows how an IR opera- 
tor O x may be built from linear combinations of direct 
products of other IR operators, O x [ and O x 2 2 , 



A Ai A 2 
A Ai A2 



)oto K2 



A 2 ' 



(Dl) 
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where A (Aj) and A (Ai) denote IR and row, 

respectively. The notation O x can refer to 

A 1 ,AZ,E\tf,t£,*'Z 1 ',* , Z 1 ",*?', or *f*. See 
the comments in the paragraph above Eq. (41) for our 
phase convention for the coefficients. 



A complete set of Clebsch-Gordan coefficients for the 
octahedral group using the basis vectors of Tables II and 
IV is given here and in Ref. [38]. The version of this 
paper that was submitted for publication was shortened 
by omitting all but a selected subset of Clebsch-Gordan 
coefficients. 

In each Clebsch-Gordan table, the resultant IR appear- 
ing on the left side of Eq. (Dl) is listed in the top row, 
and the two IRs appearing on the right side of Eq. (Dl) 
are listed in the left column. Table XVII explains how 
to read the coefficients in the Clebsch-Gordan tables in 
this appendix. 



TABLE XVII: Description for tables of Clebsch-Gordan coef- 
ficients. Squares of coefficients are listed together with their 
overall sign. 



e> Al ®e> A2 






sgn 


~C 


( A Ai A 2 \ 
VA Ai \ 2 ) 






C 


( A Ai A 2 \ 
\X Ai X 2 J 




2 



E X T 2 
E x Tl 
E 2 T{ 
E 2 T 2 
E 2 T 3 



rp\ m2 ^3 rp3 



1/4 3/4 

-1 

1/4 -3/4 

3/4 1/4 

-1 

3/4 -1/4 







G\ 


Gl H 1 


H 2 H 3 


H 4 


mGi 


H 1 


H 2 H 3 H 4 


nc\ 





1 








E 1 G\ 





-1 





T\Gl 


2/3 


1/3 





E 1 G\ 








1 


T 2 G\ 


-1/3 


2/3 





E 2 G\ 








-1 


T 2 G 2 





1/3 


2/3 





E 2 G\ 


1 








T 3 G\ 





-2/3 


1/3 







T 3 G 2 











1 


e®h 


G\ 


G\ 


(_r 2 Gr 2 n 


H 2 


H 3 H 4 






E l H l 








-1/2 1/2 












E Y B 2 


1/2 











-1/2 









E X B 3 





-1/2 








- 


1/2 






E l H 4 








1/2 








1/2 






E 2 H X 





-1/2 











1/2 






E 2 H 2 








-1/2 








1/2 






E 2 H 3 








1/2 1/2 












E 2 H 4 


1/2 











1/2 










A 2 ®E 


E 1 E 2 


A2E 1 
A 2 E 2 


-1 

1 





ml 

T$ 


m9 

T2 2 


m"3 

T 3 




A 2 ®T 2 


ml 


m9 

T 2 


m^ 

T 3 


A 2 T\ 








-1 




A2TI 








-1 


A 2 T 2 





-1 







A2T 2 





1 





A 2 T 3 


1 










A2T 3 


1 









A 2 ®Gi 


r<i r<2 

{ - r 2 U 2 


A 2 ®G 2 


G\ Gl 


A 2 G\ 
A 2 Gl 


1 

1 


A 2 G\ 
A 2 G\ 


1 
1 



A 2 ®H H 1 H 2 H 3 H A 
0-10 
1 
10 
0-100 



A 2 H l 
A 2 H 2 
A0H 3 



E®E 



E 1 E 1 
E 1 E 2 
E 2 E 1 
E 2 E 2 



E X T 2 
E'T 2 
E l T$ 
E 2 T 2 l 
E 2 T 2 
E 2 T 2 



A x A 2 E 1 E< 



1/2 1/2 
1/2 -1/2 
-1/2 -1/2 

1/2 -1/2 



E®G 2 



E 1 G 2 
E 2 G\ 
E 2 G\ 



H 1 H 2 H 3 H 4 
-1 
10 
10 
0-10 



T i T 2 T 3 T i T 2 T s 



3/4 1/4 

-1 

-3/4 1/4 

-1/4 -3/4 

-1 

1/4 -3/4 





A\ E 1 


E 2 


Tl 


T 2 


T 3 


T 1 


T| T 3 






TlTl 


1/2 














1/2 






T 1 T 2 








1/2 





1/2 









T 1 T 3 


1/3 1/6 








1/2 















T 2 T 1 





- 


1/2 





1/2 









T 2 T 2 


-1/3 2/3 
























T 2rp3 














1/2 





1/2 






T 3 T 1 


1/3 1/6 








-1/2 















T 3 T 2 














-1/2 





1/2 






T 3 T 3 


1/2 














-1/2 






m>H 


s~il r^'2 

1^1 (_ri (_r 2 


Gl 


H 1 


H 2 


H z 


H 4 


H 1 H 2 


H 3 


H 4 


tIh 1 


1/6 























-5/6 


tIh 2 


0-1/2-2/5 








01/10 








tIh 3 


1/6 





0-8/15 








0-3/10 








t\h a 


1/2 











-2/5 





01/10 





T1H 1 


0-1/3 3/5 








01/15 










1/3 





1/15 








-3/5 








T 2 R Z 


0-1/3 








0-1/15 





3/5 





TlH 4 


0-1/3 











0-3/5 





0-1/15 


TfH 1 


1/2 








2/5 








0-1/10 








TfH 2 


1/6 








8/15 





03/10 





T 3 H s 


0-1/2 











2/5 





0-1/10 


T'(H 4 


1/6 











5/6 
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7i®T 2 


A 2 


e 1 e 2 n 


n n 


rpl 
1 2 


n 


n 


nn 


1/3 


1/6 











-1/2 





1 1 1 2 








1/2 








1/2 


nn 





1/2 


-1/2 














1 1 1 2 





1/2 





1/2 








nn 


1/3 


-2/3 

















nn 








1/2 








-1/2 


nn 





-1/2 


-1/2 














nn 





-1/2 





1/2 










-1/3 


-1/6 








u 


-1/2 





Ti®G 2 


Gl 


Gl H 1 H 2 H 3 H 4 




















-1 














2/3 





1/3 












TfG* 


-1/3 





2/3 












7?G| 





1/3 2/3 















TfG£ 





-2/3 1/3 




















0-1 















T 2 <g)r 2 


Ax 


E 1 E 2 n 


n n 


rpl 

1 2 


n 


n 


nn 





0-1/2 











-1/2 





nn 








-1/2 








-1/2 


nn 


1/3 


1/6 


-1/2 














nn 








1/2 








-1/2 


nn 


1/3- 


2/3 

















nn 





1/2 





1/2 








nn 


1/3 


1/6 


1/2 














nn 





-1/2 





1/2 








nn 





-1/2 











1/2 





T 2 ®Gi 


Gl 


G\ H 1 H 2 


H 3 H 


4 










T^GJ 





-2/3 1/3 
















T\G\ 





0-1 
















nc\ 


1/3 





-2/3 










na\ 





-1/3 -2/3 
















T£G\ 








1 













nc\ 


-2/3 





-1/3 










T 2 ®G 2 


G\ 


G\ H 1 H 2 


H 3 H 4 












T\G\ 





-2/3 - 


1/3 










T\G\ 








-1 










nc\ 


1/3 


2/3 













not 





-1/3 


2/3 










nc\ 





1 













nc\ 


-2/3 


1/3 














Tjg>H 


G\ 


Gl 


Gl 


Gl 


H 1 


H 2 


H 3 H 4 H 1 


H 2 


H 3 H 4 


nn x 








1/2 

















1/2 


T\n 2 











1/6 


2/3 





01/6 








nH 3 


1/2 

















0- 


1/2 





nn 4 


0- 


1/6 











0- 


2/3 





1/6 


T 2 H X 


0- 


1/3 














1/3 





1/3 


T 2 H 2 








1/3 











01/3 





0-1/3 


T 2 H 3 











1/3- 


1/3 





01/3 








T 2 H A 


-1/3 











0- 


1/3 


0- 


1/3 





T^H 1 


1/6 











0- 


2/3 





1/6 





T 3 H 2 


0- 


1/2 

















0- 


1/2 


TiH 3 





0- 


1/6 











02/3 





1/6 


T 3 H 4 








0- 


1/2 








01/2 









G\G\ 
G\G\ 
G 2 G\ 
G\G\ 



G 2 (8)G 2 



2^2 

G 2 G 2 
r<2rii 

{ ~ r 2 U 2 

r<2ri2 



10 

1/2 1/2 

•1/2 1/2 

1 



Ax n n n 



10 

1/2 1/2 

-1/2 1/2 

1 



Gi®G 2 


A 2 


rpl 
1 2 


T 2 
1 2 


n 


G\Gl 











-i 


Lr l Lr 2 


1/2 





-1/2 





r<2r<~t 

Lt 1 Lt 2 


-1/2 





-1/2 





G\G\ 





1 











E 1 


E 2 


n 


n 


n n 


n n 


GlH 1 


1/2 











1/2 


G\H 2 








1/4 





3/4 





G\H 3 


1/2 








1/2 








G\H 4 














3/4 


1/4 


GlH 1 








-3/4 





1/4 





G 2 H 2 


1/2 








-1/2 








G 2 H 3 














-1/4 


3/4 


G 2 H 4 


1/2 











-1/2 



G-2&H 


E 1 


e 2 n 


n 


If 


rpl 
1 2 


T 2 

1 2 


n 


GlH 1 


1/2 





1/2 














GlH 2 











-3/4 








-1/4 


GlH 3 





-1/2 











-1/2 





GlH 4 





1/4 








3/4 








GlH 1 











-1/4 








3/4 


GlH 2 





-1/2 











1/2 





GlH 3 





3/4 








-1/4 








G 2 2 H 4 


1/2 





-1/2 

















A x 


A 2 


E 1 


E 2 


n 


T 2 


Tf T 1 


n 


nn 


T 2 


nn 


T 2 


n 


H'H 1 


























25 








3 











H l H 2 





1 





1 























-1 








1 





H X H 3 














3 








3 








1 








1 








H 1 !! 4 


1 





1 








9 








-1 























H 2 H X 





1 





-1 























-1 








-1 





H 2 H 2 














-4 








9 








3 

















H 2 H 3 


-1 





1 








-1 








-9 























H 2 H 4 




















3 








3 








-1 








1 


H^H 1 














3 








3 








1 








-1 








H 3 H 2 


1 





-1 








-1 








-9 























H 3 H 3 




















-4 








9 








-3 











H 3 H 4 





-1 





1 























-1 








-1 





H 4 H X 


-1 





-1 








9 








-1 























H 4 H 2 




















3 








3 








-1 








-1 


H 4 H 3 





-1 





-1 























-1 








1 





H 4 H 4 




















25 








-3 

















denom 


4 


4 


4 


4 10 20 10 40 20 40 


8 


4 


8 


2 


4 


2 



Numerators given in a column of this table should be divided 
the listed denominator before taking the square root. 



